Propagation of Coulomb-correlated electron-hole pairs in semiconductors with 

correlated and anticorrelated disorder 
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Local ultrafast optical excitation of electron-hole pairs in disordered semiconductors provides 
the possibility to observe experimentally interaction-assisted propagation of correlated quantum 
particles in a disordered environment. In addition to the interaction driven derealization known 
for the conventional single-band TIP-(two-interacting-particles)-problem the semiconductor model 
has a richer variety of physical parameters that give rise to new features in the temporal dynamics. 
These include different masses, correlated vs. anticorrelated disorder for the two particles, and 
dependence on spectral position of excitation pulse. 

PACS numbers: 72.15.Rn, 73.20.Dx, 73.20.Jc, 78.47.+P, 78.66.-w 



Introduction. It is by now well established that the 
center-of-mass (COM) motion of two interacting par- 
ticles (TIP) in a single tight-binding band of a one- 
dimensional Hamiltonian possessing diagonal disorder is 
spatially more extended than both their relative motion 
and the motion of the particles without interaction |Q. 
Although this model is a fundamental paradigm for the 
more general question of the interplay of disorder and 
interactions in many-body systems, it suffers from not 
being susceptible to experimental verification. Further- 
more, it is not clear whether this finding, relevant for two 
particles, can be generalized to real transport situations 
in dense Fermi systems. 

On the other hand, low-intensity optical excitation 
in disordered semiconductor systems produces strongly 
Coulomb-correlated electron-hole pairs with large mutual 
separation such that their motion can be considered to 
be that of an essentially isolated pair. Optical pump- 
probe experiments are in principle suited to record the 
coherent dynamics of the pair after pulsed excitation on 
a ps-time scale before scattering with acoustic phonons 
becomes relevant. 

Following this idea the temporal traces of the partic- 
ipation number of the electron and the hole have been 
studied in a previous paper ^| by solving the equation 
of motion for the two-particle amplitude Pij(t) = (djCj), 
which is the interband coherence related to the optical 
polarization. Here di (dj) are hole (electron) operators 
at site i (j). In the low intensity (with respect to the 
exciting light pulse) limit the sum rules n\j = X]; PijPu 
and riy = J2 t PjiP*i yield the electron (hole) intraband 
quantities. Their diagonal elements n^ h are the time 
dependent densities of electron and hole. Defining the 
participation number A(t) = (X^ n ri) _1 the evolution of 
the excited electron and hole wave packet can be plotted 
as a function of time. Alternatively, the COM-coordinate 
R and the relative coordinate p of the electron-hole wave 
packet can be calculated directly from the pair ampli- 



tudes. 

It has been demonstrated that the interaction induces 
an enormous change of the dynamics if compared to the 
non-interacting case. Instead of an exponential rise of 
A(t) towards a saturation value on the time scale of the 
excitation pulse, in the interacting case a slow, diffusion- 
like rise of A is seen that does not seem to saturate in the 
limited time regime accessible to the numerical calcula- 
tion (and relevant for the real physical situation). For 
brevity this feature will be called "enhancement" in the 
following in accordance with the notion of other work in 
the field, although it is not implied that from the present 
dynamical calculation in a finite time domain anything 
like the enhancement of a localization length can be de- 
duced. The dependence of the enhancement on the sign 
of the interaction has been studied and it was found that 
the dependence on the sign vanishes only for excitation 
in the center of the continuum. 

In this contribution the effect of correlated versus an- 
ticorrelated disorder and of the interaction strength is 
studied. 

Model. A two-band Hamiltonian Hq is considered 
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with nearest neighbour coupling J e (J h ) for electrons 
(holes) and diagonal disorder given by a box-shaped dis- 
tribution function of the single-site energies ef (e^ 1 ) of 
total width W . The N sites form a linear chain with 
spacing a and periodic boundary conditions. 

The interaction is given by a regularized Coulomb po- 
tential in monopole-monopole form 
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with 

U e 2 

Va = ~ (3) 

Ait££q m + a 

where a = 5a and U is a dimensionless parameter quan- 
tifying both strength and sign of the interaction. 

The optical excitation is given by a local dipolar cou- 
pling to the light field 

N 

Hi = -^{inE*{t)diCi + h.c). (4) 
i=l 

For local excitation at site the optical dipole matrix 
elements are taken to be [ii = /ito^i,o- The light field is 
given by 

E(t) = {Tt)- 1/2 a- 1 exp (-(t/a) 2 ) exp (-iu Q t) (5) 

with central frequency luq and temporal width a. 

Equation of Motion. The equation of motion for the 
pair amplitude pij is derived using the total Hamiltonian 
H = Hq + He + Hi . In the low intensity limit linear re- 
sponse theory is valid and the equation of motion reads 
(with h = 1) 

N 

+ i^2{J e pu + J h pij) + injEtySij. (6) 
i=i 

This equation of motion is solved numerically for M re- 
alizations (typically M — 20. ..40) of the disorder drawn 
from the distribution of site energies. The observables are 
then configurationally averaged over these realizations. 

Choice of Parameters. As the present model de- 
scribes a situation that is in principle accessible to ex- 
periments, the interaction strength is chosen to yield 
exciton binding energies in the range of values typi- 
cal for low-dimensional semiconductor heterostructures. 
The intraband couplings J e ' h and the strength of dis- 
order W e ' h are free parameters. However, the tight- 
binding model is thought to model the conduction and 
valence band extremities of a disordered semiconductor 
heterostructure. In this communication we show data 
for J e = J h = 20meV exclusively, i.e. equal electron and 
hole masses are assumed. The lattice constant, equal to 
a disorder length scale, is taken to be a = 20A. The 
external optical pulse resembles typical laser pulses used 
in ultra-fast optical experiments on coherent phenomena. 
Its central frequency luq is chosen to be situated below 
the center of the absorption band calculated without in- 
teraction, in order to model the dynamics of electron- 
hole pairs in the continuum in a range of energies that 
excludes LO-phonon emission. Specifically, the gap of 
the noninteracting case is the origin of the energy axis, 
the band center is at A = 80 meV, and ojq = 40 meV. 



For this excitation condition the sign of the interaction 
is not irrelevant |S|. Here we use attractive interaction 
throughout. A temporal pulse width a =100 fs is used 
corresponding to a spectral width of 22 meV. The num- 
ber of sites N is taken large enough (typically N = 240) 
such that for not too small disorder the locally excited 
wave packet does not reach the boundary of the sample in 
typically some ten picoseconds. After this time acoustic 
phonon scattering will lead to effective dephasing and the 
coherent phenomena studied here will be destroyed. For 
principal studies, however, the calculation is sometimes 
performed for much longer times, although the results 
are no longer relevant to experiments. 

Observables. The extension of the electron and hole 
wave packets as a function of time can be character- 
ized by the above mentioned participation number A e ' h 
which arc single-particle quantities. Using the pair am- 
plitude pij, the two-particle quantities COM-coordinatc 
R = (J2 i ■ \pij \ 2 (i + j) 2 /2) 1 / 2 and relative coordinate 

P = j \Pi]\ 2 (i ~ i) 2 /2)^ 2 can be defined. Their ratio 
R/ p is a measure of the interaction induced enhancement. 

Correlated and Anticorrelated Disorder. While 
previous results [Q] have been obtained using uncorre- 
cted disorder, here the influence of correlation is stud- 
ied. Disorder for electrons and holes is called correlated 
if the site energies of electrons equal that of the mirror 
image of the hole energies with respect to the gap cen- 
ter. This kind of disorder is expected in heterostructures 
due to interface roughness leading to a spatially fluctu- 
ating confinement of single-particle wave functions. On 
the other hand, if the energy separation within an iso- 
lated site remains constant, i.e. equal to A, the disor- 
der is called anticorrelated. A spatially fluctuating elec- 
trostatic field would induce this kind of disorder. The 
present model resembles that of the conventional single- 
band-TIP if electrons and holes have equal masses and 
if correlated disorder (with W e = W h ) is applied. Elec- 
tron and hole wavefunctions are then pairwise equal, both 
particles live in the same environment. For U — all op- 
tical matrix elements connecting these pairs of states are 
equal, all others are zero. Therefore, the optical spec- 
trum for U = resembles the density of single-particle 
states, however, with bandwidth being the sum of that 
of the two bands. As in the last years the conventional 
TIP-model has been widely studied we here concentrate 
mainly on the anticorrelated situation. 

Figure [j] shows COM-traces for correlated and anti- 
correlated disorder for U = and U = 3. Remark- 
ably, the enhancement for anticorrelated disorder is much 
more pronounced compared to the correlated case. These 
traces have been calculated for an optical excitation en- 
ergy ujo = 40 meV. Looking at the optical spectra for 
the anticorrelated case, Fig. |^, one realizes a peak that 
for U = Q lies in the center of the band at A = 80 meV 
and shifts towards lower energy for increasing (attrac- 
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tive) U. For the excitation condition leading to the large 
COM-enhancement at U = 3 pairs of states in the vicin- 
ity of this peak are excited. The origin of the peak can 
be traced back to optical transitions connecting states in 
the tails of the single-particle bands. Even for moderate 
disorder these pairs of states with the nearly identical 
transition energy are strongly localized at the same po- 
sition, and therefore their optical matrix elements are 
large. Exact eigenvectors for a short sample (N = 10) 
have been calculated for U = and U = 3, confirming the 
strong overlap within the contributing pairs. Obviously, 
these strongly localized tail states are not responsible for 
the large enhancement. 

The fact that the enhancement for anticorrelated dis- 
order is larger than that for correlated disorder is not 
fully understood yet. A prominent difference exists in 
the equations of motion for py for the two cases. For the 
anticorrelated case the disorder is absent in the equation 
for the diagonal elements pa, while it is present in the 
respective equation for the correlated case. 

The behavior of excitons in disordered systems is often 
H discussed in terms of the relative coordinate r = i — j 
and the COM-coordinate x = (i + j) /2 instead of the in- 
dices i and j (for equal masses, not to be confused with R 
and p). After the transformation — > (r, x) the rela- 
tive coordinate r can be integrated out in the equation of 
motion if, e.g., the disorder W is smaller than the exciton 
binding energy and if discrete excitonic resonances are 
considered. In our case the equation of motion in terms 
of x would not contain any disorder in the anticorrelated 
case, in contrast to the correlated case. Although for the 
parameters used in Fig. [I] this approach is not strictly 
valid, it suggests a possible solution for the problem at 
hand. 

Finite-Time Scaling. In order to quantify the en- 
hancement finite-time scaling is applied in the following 
sense. Although with a dynamical calculation for finite 
times it is impossible to decide whether the two-particle 
packet is localized or not a long-time saturation value of 
the COM-coordinate Roo is assumed for practical pur- 
poses. While the temporal rise of the interaction-free 
traces is exponential with time scale given by that of the 
pulse duration, the short time behavior in the interact- 
ing case looks like a diffusive process. So the following 
interpolation relation is used: 



R(t) = {(Dt) 



-1/2 



(7) 



Figure |^ shows the very good quality of typical fits for 
times larger than the pulse width. It is found that while 
the diffusivity D is weakly dependent on disorder and in- 
teraction, Roo shows a much stronger dependence. In the 
following Roo is taken to characterize the enhancement. 

Contour Plots. Figure ^ shows contour plots of the 
two-particle wave function, |py| 2 , for anticorrelated and 
correlated disorder at time t = 165 ps. All other pa- 



rameters are equal to that used in Fig. 13. Again, the 
anticorrelated situation shows much larger enhancement. 

Dependence of Enhancement on Interaction 
Strength. Recently the dependence of the enhancement 
on interaction U has been discussed and a duality rela- 
tion describing the behavior for large and small U has 
been proposed |5[. The present calculations yield similar 
results, see Fig. | for anticorrelated disorder. The ratio 
Roo/ Poo is close to unity for small U , while it has a max- 
imum around U = 3. ..4 that increases with increasing 
disorder. Of course, both Roo and poo decrease with in- 
creasing disorder, however, poo faster than Roo- Remark- 
ably, the position of the maximum does not depend on 
disorder. This behavior of Roo/ Poo is apparently related 
to the fact that for the chosen excitation energy loq the 
center of the optical spectrum, indicated by the above 
mentioned peak, shifts (c.f. Fig. 2) through the spec- 
tral position of the laser pulse with changing U. Around 
U = 3. ..4 the peak coincides with the excitation energy 
u)o, while at larger U the excitation takes place again in 
a region where states outside the center of the single- 
particle bands are optically coupled. It is known that 
the enhancement is more pronounced for less localized 
single-particle states QJ^]. These states are situated in 
the center of the single-particle bands. The Coulomb en- 
hancement is therefore largest if the electron-hole packet 
is generated in the center of the optical spectrum since 
the interaction couples these weakly localized states most 
effectively. 

Conclusions. For a situation that is in principle ac- 
cessible to ultra-fast coherent optical experiments, the 
delocalizing action of the Coulomb interaction on opti- 
cally generated electron-hole pairs has been studied. The 
temporal traces of the center-of-mass coordinate R and 
the relative coordinate p have been fitted by an interpola- 
tion scheme describing diffusive and localized behavior at 
small and large times, respectively. It should be noted, 
however, that the calculation presented here is unable 
to yield definite answers about localization or otherwise. 
Correlated and anticorrelated disorder is studied and it 
is found that the enhancement of Roo/ Poo is much more 
pronounced for anticorrelated disorder. 

Acknowledgements This work is supported by DFG, 
SFB 383, the Leibniz Prize, OTKA (T029813, T024136, 
F024135), SNSF (2000-52183.97), and the A. v. Hum- 
boldt Foundation. 



* Permanent address: Department of Theoretical Physics, 
Institute of Physics, Technical University of Budapest, H- 
1520 Budapest, Hungary 

[1] D.L. Shepelyansky, Phys. Rev. Lett. 73, 2607 (1994). 



3 



[2] D. Brinkmann, J.E. Golub, S.W. Koch, P. Thomas, K. 

Maschke, and I. Varga, Eur. Phys. J.B. 10, 145 (1999) 
[3] L. Banyai, I. Galbraith, C. Ell, H. Haug, Phys. Rev. B36, 

6099 (1987). 

[4] S.D. Baranovskii and A.L. Efros, Fiz. Tekh. Poluprovodh. 

12, 2233 (1978) [Sov. Phys. Semicond. 12, 1328 (1978)] 
[5] X. Waintal, D. Weinmann, and J.-L. Pichard, Eur. Phys. 

J.B. 7, 451 (1999) 

FIG. 1. Center-of-mass (R) traces for U = and 

17 = 3 for anticorrelated and correlated disorder. Disorder 
W e = W h = w = 6Q me y^ 



FIG. 3. Fit of the center-of-mass trace by Eqn. (7) for 
different interaction strengths U. Other parameters are the 
same as in Fig. 1. 



FIG. 4. Contour plots of \pij\ 2 at 165 ps for anticorrelated 
(left) and correlated (right) disorder. [7 = 3. Other parame- 
ters are the same as in Fig. 1. 



FIG. 2. Absorption spectra for U = and 17 = 3 for cor- 
related (dotted line) and anticorrelated (solid line) disorder. 
Other parameters as in Fig. 1. 



FIG. 5. Dependence of the enhancement Rao/poo on in- 
teraction strength U for three different disorder strengths 
W e = W h = W. Inset shows R x (full symbols) and p^ 
(open symbols) vs. U separately. 
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